Equation of state of hypernuclear matter: impact of hyperon— scalar-meson couplings 
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We study the equation of state and composition of hypernuclear matter within a relativistic density 
functional theory with density-dependent couplings. The parameter space of hyperon-scalar-meson 
couplings is explored by allowing for mixing and breaking of SU(6) symmetry, while keeping the 
nucleonic coupling constants fixed. The subset of equations of state, which corresponds to small 
values of hyperon-scalar-meson couplings allows for massive M < 2.25Mq compact stars; the radii 
of hypernuclear stars are within the range 12-14 km. We also study the equation of state of hot 
neutrino-rich and neutrinoless hypernuclear matter and confirm that neutrinos stiffen the equation 
of state and dramatically change the composition of matter by keeping the fractions of charged 
leptons nearly independent of the density prior to the onset of neutrino transparency. We provide 
piecewise-polytropic fits to six representative equations of state of hypernuclear matter, which are 
suitable for applications in numerical astrophysics. 
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I. INTRODUCTION 

The integral parameters of compact stars depend on 
their equation of state (hereafter EoS) at high densities. 
Measurements of pulsar masses in binaries provide the 
most valuable information on the underlying EoS because 
these, being deduced from binary system parameters, are 
model independent within a given theory of gravity 
The recent discovery of a compact star with a mass of 
1.97 Mq (PSR J1614-2230) which sets an observationally 
"clean" lower bound on the maximum mass of a compact 
star via the measured Shapiro delay spurred an in- 
tensive discussion of the phase structure of dense matter 
in compact stars consistent with this limit [3l-l27j . Pulsar 
radii have been extracted, e.g., from modeling the X-ray 
binaries under certain reasonable model assumptions, but 
the uncertainties are large [28l - [3]| . An example, that we 
will use in our analysis, is the pulse phase-resolved X- 
ray spectroscopy of PSR J0437-4715, which sets a lower 
limit on the radius of a 1.76 Mp solar mass compact star 
i? > 11 km within 3cr error |31| . 

Large masses and radii are an evidence for the rela- 
tive stiffness of the EoS of dense matter at high densi- 
ties. Large densities may require substantial population 
of heavy baryons (hyperons), because these become en- 
ergetically favorable once the Fermi energy of neutrons 
becomes of the order of their rest mass. Their onset then 
reduces the degeneracy pressure of a cold thermodynamic 
ensamble and softens effectively the EoS of dense mat- 
ter. This decreases the maximum mass of a compact 
stars to values which contradict the observation of mas- 
sive compact stars in Nature. This controversy between 
the theory and observations is the essence of the "hyper- 
onization puzzle" in compact stars. 

Although the emergence of heavy baryons (mainly 
and A hyperons) was considered even before the discovery 
of pulsars and their identification with neutron stars js^ , 
their existence in the cores of neutron stars is still elu- 
sive. Seminal work on interacting hyperonic matter was 
carried early after the discovery of pulsars [1^ [s^] . Sys- 



tematic studies of interacting dense hypernuclear matter 
were carried out with the advent of relativistic density 
functional methods [sBl - ii^ l , for reviews see [1^, . Most 
of the these works predicted masses that are not much 
larger than the canonical mass of a neutron star (in con- 
tradiction with modern observations) . Masses on the or- 
der of < 1.8Mq were obtained later on in non- relativistic 
phenomenological models [l^ . Fully microscopic models 
based on hyperon-nucleon potentials, which include the 
repulsive three-body forces, predict low maximal masses 
for hypernuclear stars [i^, . 

Finite temperature hyperonic matter in the presence 
of thermal bath of neutrinos may differ significantly from 
cold hypernuclear matter found in evolved neutron stars. 
Previous studies [50l - l53| showed that the EoS becomes 
stiffer, the deleptonizing effect of charged hyperons at 
zero temperature is removed in the presence of neutri- 
nos and the sequence in which charged hyperons appear 
changes. 

In this work we contribute to the investigation of "hy- 
peronization puzzle" by studying the parameter space 
of the couplings of hyperons to the light scalar octet of 
mesons. The parameter space is partially motivated by 
the hypernuclear potentials as well as investigations of 
the external field QCD sum rules that involve SU(6) sym- 
metry breaking and mixing (s^ . [ssj . Because the struc- 
ture and even content of the scalar mesons is still a sub- 
ject of discussion it appears necessary to investigate this 
particular domain of hypernuclear physics with respect 
to the consistency with the compact star observations. 
Conversely, as we show below, the considered parameter 
space can in fact account for massive hyperonic compact 
stars, thus providing a possible solution to the "hyper- 
onization puzzle" . 

To study the EoS of hypernuclear matter we use a 
model based on relativistic density functional theory of 
Walecka and a density dependent parameterization of 
the nucleon-meson couplings of Rcf. (56|] . The EoS and 
composition of matter will be studied both in the zero 
temperature limit relevant to mature compact stars, as 
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well as at finite temperatures relevant for their hot proto- 
neutron star stage of evolution. In the latter case we as- 
sume neutrino-rich matter in /3-equilibrium and explore 
the effects of finite temperature and neutrino content on 
the stiffness and composition of matter. We confirm the 
general features found in the previous analysis fsol - fssj 
and provide a detailed account of the EoS and composi- 
tion of matter within our model. 

This paper is structured as follows. In Sec. [ll] we dis- 
cuss the theoretical set-up of the relativistic density func- 
tional theory. The parameterization of the coupling con- 
stants of the theory is discussed in Sec. IIIII Sec. |IV] is 
devoted to the discussion of the choice of the hyperon- 
meson coupling constants. We present our results in 
Sec. |Vl Finally, our conclusions are collected in Sec. IVD 



II. THEORETICAL MODEL 



stellar matter and form a statistical ensamble in equi- 
librium. In mature compact stars at zero temperature 
the neutrinos do not contribute to the thermodynamical 
quantities of matter because they are out of equilibrium 
and their chemical potential vanishes. 

The Lagrangian density (P) yields the pressure and the 
energy density, which at finite temperature read respec- 
tively 
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The relevant degrees of freedom in nuclear matter at 
high density are nucleons and hyperons. In particular, 
we choose as degrees of freedom of our relativistic mean 
field model the baryon octet (J^ = i^), and we study 
the interaction of these fields with the isoscalar-scalar, 
tr, isoscalar-vector, w^, and isovector-vcctor, p^, mesons. 
Besides these fields, we also consider the presence of lep- 
tons, e~ , fi~ . In the case of finite temperature matter, 
which describes proto-neutron star matter, the model has 
the neutrino degrees of freedom, because they are ex- 
pected to be trapped in a proto-neutron star within the 
first minute after its birth. The introduction of additional 
variables, the neutrino chemical potentials, requires addi- 
tional constraints, which we provide by fixing the lepton 
fractions, Yu, appropriate for conditions prevailing dur- 
ing the evolution of a proto-neutron star [5l| . 

The relativistic Lagrangian reads 
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where the _B-sum is over the = baryon octet, ipB 
are the baryonic Dirac fields with masses niB- The in- 
teraction is mediated by the cr, and meson fields 
with the uj^i, and p^i, field strength tensors and masses 
TOcr, nii^, and mp. The baryon- meson coupling constants 
are denoted by gBBm (their numerical value will be dis- 
cussed in Sec. |llT] below) . The A-sum in Eq. ([T]) runs 
over the leptons e~,/x~,i/e and i/^ with masses mx and 
the last term is the electromagnetic energy density. The 
contribution of neutrinos to the sum above is included 
only at non-zero temperature, when they are trapped in 
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where cr, cjq and po3 sxc the nonvanishing mesonic mean 
fields, Jb is the baryon degeneracy factor, m*^ = ms — 
gaBB is the effective baryon mass, /i^ = ^b — guBB^o — 
gpBBhPz is the baryon chemical potential including the 
time component of the fermion self-energy, 1^ is the third 
component of baryon isospin and /(y) = [l-|-exp(?//T)]~"'^ 
is the Fermi distribution function with T being the tem- 
perature. 

At non-zero temperature the net entropy of the matter 
is the sum of the baryon contribution 
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and the lepton contribution 
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where f{y) = 1 — f{y). The free-energy density is then 
given by 



F = e-T{SB + SL). 



(6) 
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III. DENSITY-DEPENDENT 
PARAMETRIZATION DD-ME2 

Below, we shall work with a density dependent pa- 
rameterization, which is designed to account in an eco- 
nomical manner the many-body correlations that arise 
beyond the mean-field approximation. The density de- 
pendence of the coupling constant accounts for the in- 
fluence of the medium on the scattering of baryons. We 
choose to work with the DD-ME2 parameterization, de- 
signed in Ref. [56j . which introduces an explicit density 
dependence in the nucleon-meson couplings gNNa, 9nnui 
and gNNp- The phenomenological Ansatz for the density 
dependence for the u and meson coupling constants 
is 



gNNiips) = gNNi {po)hi{x), i 
where pB is the baryon density. 



(T, OJ, 



hi{x) 



'l + Ci{x + diy 



(7) 



(8) 



is a function oi x ~ pb/po, Po is the nuclear satura- 
tion density. The parameters in ([5]) are not independent. 
The five constraints hi{l) = = hZ{l),h'/{0) = 0, 

reduce the number of independent parameters to three. 
Three additional parameters in the isoscalar-scalar chan- 
nel are gNNa{po), gNNu{po) and the mass of the 
phenomenological a meson. Microscopic calculations 
show that the p^-meson coupling vanishes at high den- 
sities [131 • Therefore, instead of the Ansatz (|8]) one can 
use a simpler exponential form of density dependence 



gNNpips) = gNNpipo) exp[-ap(x - 1)]. 



(9) 



Thus, the parameterization we will be using has in total 
eight parameters, which are adjusted to reproduce the 
properties of symmetric and asymmetric nuclear matter, 
binding energies, charge radii, and neutron radii of spher- 
ical nuclei. The parameters of the DD-ME2 effective in- 
teraction are shown in Table HI 



TABLE I: Meson masses and their couplings to the baryons 
in DD-ME2 parametrization. 





a 


U) 


P 


rrii [MeV] 550.1238 783.0000 763.0000 


QNNiipo) 


10.5396 


13.0189 


3.6836 


ai 


1.3881 


1.3892 


0.5647 




1.0943 


0.9240 




Ci 


1.7057 


1.4620 






0.4421 


0.4775 





When the coupling constants are density dependent, 
the thermodynamical consistency (specifically the energy 
conservation and fulfillment of the Hugenholtz-van Hove 
theorem) require the inclusion of the so-called rearrange- 
ment self-energy [ssj . This contributes to the pressure. 



but not to the energy of the system. The pressure be- 
comes 

Pr = P + PB^r, (10) 

where the rearrangement self-energy, S,., is given by 



dgNNuj dgNNa 
= — WoPS CrpS: 



dpi 



dpi 



(11) 



where ps is the scalar density. It can be verified that the 
contribution from the rearrangement self-energy restores 
the thermodynamical relation 



dpB \Pb 



(12) 



which is otherwise violated. Below we shall always con- 
sider the pressure including the rearrangement term and 
will drop the subscript r. 



IV. HYPERON-MESON COUPLING 
CONSTANTS 

To extend the description of matter to full baryon octet 
we need the hyperon-meson coupling constants. Because 
the information on the properties of hypernuclear matter 
is far less extensive than for nucleons their values can not 
be fixed with certainty and several approaches have been 
used in the literature. 

One way to fix these couplings is to consider the SU(3)- 
flavor symmetric model, namely the octet model. In this 
model the degrees of freedom are represented by the low- 
est non-trivial irreducible representation (IR) of the sym- 
metry group which is physically possible, {8}, that is the 
baryon octet and the mesonic octet. 

Let us consider first the interaction between the baryon 
octet, = and the vector meson octet, = 
1^. Within the assumption of SU(3)-flavor symmetry, 
one can express the meson-baryon coupling constants in 
terms of only two parameters [59| , the nucleon-p^ meson 
coupling constant gNNp and the F/ {F + D) ratio of the 
vector octet ay: 

gNNp = 38, ffHSp = -58(1 - 2av), 

gY,T.p = S^gay, .9AAp = 0, 

gNNuf, = -^5'8(4av - 1), g^.^^^ = --^58(1 + 2av), 

2 . ^ 2 ^ ^ 

gsswg = - av), gwius = --j=g8, (1 - ay), 



(13) 

where we have shown only the coupling constants rele- 
vant for our model. The vector-meson dominance model 
[60j predicts ay = 1, which is a result of the universal 
coupling of the meson to the isospin current. This 
result is also in agreement with recent QCD sum rules 
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calculations [55|. Therefore, we set ay = 1 in hyperon- 
meson coupling constants in Eq. p^ . i.e., 



5HHp 



58, 



5SEp = 258, ffAAp = 0. 



(14) 



To describe the mixing between the singlet and the octet 
members of the vector nonet, one has to introduce two 
additional parameters - the flavor singlet coupling con- 
stant, gi, and the vector mixing angle, 9v- Therefore, 
the coupling of the baryon to the physical w^-meson state 
now reads 



gBBu, = cos9vgi + sinOv gBBu 



(15) 



It is widely accepted that the mixing between non- 
strange and strange quark wave functions in the uj and (f) 
mesons is ideal |6ll. |63| . In this limit tstnOv = l/v^ and 
assuming again ay = 1 we obtain the baryon-w^ meson 
couplings as 



gNNu) 

sin 9v 

sin 9v 



V2gi 
V2g, 

gAAu 



sin 9v sin 6v 



Vigs, 
Vigs, 
= V2gi. 



(16) 



Next we assume that the nucleons do not couple to pure 
strange mesons ((/f^^). One then obtains the relation be- 
tween the singlet and the octet coupling constants [gI] : 
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= Vg 



58, 



(17) 



which allows one to eliminate gi from Eqs. ([T5|) . Thus, 
the hyperon-oj^ meson coupling constants can be related 
to the nucleon one 



1 



5E 



:gNNu 



ghAu 



:gNNu 



(18) 



We are now left with the hyperon-scalar meson (a) cou- 
pling constants. In order to fix these parameters, we 
consider the relations among the coupling constants for 
the scalar octet [H^l , which are given [in analogy with the 
relations for the vector meson octet in Eq. (fT3|) ] by 



gNNao ~ 5S, gSEao 
g^^aa = '^gsas, 5AAao 



^9s(l - 2as), 



gNNas 
5SSO-8 



1 

7! 

2 

V3 



g 3(4.0:3 - 1), 



= 0, 

5HH<T8 



5s(l-as), 5AAo 



1 
2 

71 



.gs(H-2as), 

.gs(l - as). 

(19) 



The mixing between the singlet and the octet states of 
the scalar nonet is given by 



gBBa = COS^sgi + SmtlsgBBa 



(20) 



Then, the a meson-baryon coupling constants are 

gNNa = COS 6*5 gi + sin 6*5 (4as - l).gs/\/3, (21) 

5AA<T = cos 6*5 gi - 2 sin 6*5 (1 - Q!5)5s/V3, (22) 

5ES<T = cos 6*5 gi + 2 sin 6*5 (1 - Q!5)5s/%/3, (23) 

gssa - 008 03 51 - sin 61s (1 + 2as)gs/V3. (24) 

Combining the couplings, we obtain 

gAAcr+g^^a = 2 COS 6*5 51, (25) 

gNNa+gssa = 2 COS 03 gi - 2 sm0 3 (1 - as)5s/\/3 
= cos 6*5 gi +gAAa- (26) 

Eliminating gi from these expressions we obtain a rela- 
tion 



'2igNNa + gSEa) = S^AAo- + 5SSo-, 



(27) 



which is valid for arbitrary values of the four parameters 
cks, 51, 9s and ^5. In particular, it is satisfied for the 
values of the coupling constants in the SU(6) symmet- 
ric quark model. Indeed, the latter model assumes that 
the a meson is a pure up and down state, therefore for 
strangeness —1 A and E baryons^AAo- = 5sS(t ~ ^gNNa, 
whereas for strangeness — 2 S baryons 55=0- = \gNNG- 
Due to the constraint (P7)) the parameter space of our 
model is spanned by three out of the four parameters 
Q^s, 5i, 9s and 03. Further constraints on the param- 
eter space can be placed on account that the hyperon 
coupling constants must be positive and less than the 
nucleon coupling constant. This implies two additional 
constraints, which can be translated into constraints on 
the range of variability of the hyperon-scalar meson cou- 
pling constants. Indeed, by expressing one of the hyperon 
coupling constants, say gsscn in terms of the others with 
the help of Eq. (gT]) one finds 

g-EEa = ^{igAAa + gY.y.a) " 9NNa, (28) 



and therefore, by requiring that 

< 5hH(t < 9NNa, 



(29) 



one is left with the following two simultaneous inequali- 
ties 



9NNa < -j^iigAAa + .gsSo-) < ^g^No 



(30) 



To explore the parameter space spanned by the hyperon- 
a meson couplings we use the following approach. We 
fix the value of ^aAo- at the value provided by the Ni- 
jmegen Soft Core (NSC) hypernuclcar potential model of 
Ref. [131 and vary the range of couplings gsT.a within the 
limits provided by Eq. (pO)). The corresponding gsEa 
couplings are found from Eq. (pS)) . Analogously, to ex- 
plore the range of admissible ^aAo- couplings, we next fix 
the value of gsso- to the value provided by the NSC model 



and change the gp^Ka couphng within a range which keeps 
the values of g^B.^ consistent with Eq. (P^ . The param- 
eter space is further limited by the requirement that the 
maximal values of _gAA(T and gssa should be below their 
values in the SU(6) symmetric model. This latter con- 
straint is motivated physically by our search for a stiff 
EoS of hypernuclear matter. 



V. RESULTS 

The matter in evolved compact stars is charge neu- 
tral and in /3-equilibrium. These two conditions are im- 
posed when computing the EoS at zero temperature. For 
young neutron stars at non-zero temperature a thermal 
ensamble of neutrinos should be included for tempera- 
tures T > 5 MeV. In this case a common practice is to 
vary the lepton fraction of matter in a certain range com- 
patible with astrophysical simulations. We shall treat 
below these distinct cases in turn. 



A. Zero temperature 

We start the discussion of our results by comparing the 
EoS of purely nucleonic matter, computed using density- 
dependent (DD) ME2 parameterization of nuclear mat- 
ter, to the EoS where interacting hyperons are added 
with couplings taken according to SU(6) quark model, 
see Fig. [TJ The appearance of hyperons, which is trig- 
gered by the fact that the cost of having a hyperon heav- 
ier than a nucleon is energetically more favorable than a 
neutron at the top of the Fermi sea, always softens the 
EoS. Because in our search the coupling constants are 
bound from above by their values in the SU(6) model, 
these two EoS correspond to the stiffest and softest EoS 
in our collection. 

Our next step is to explore the impact of the varia- 
tion of the hyperon-scalar meson coupling on the EoS. 
As explained in the previous Sec. IIV[ we will keep the 
hyperon-vector meson couplings fixed at the values given 
by Eqs. ([H]) and and vary the hyperon-scalar me- 
son coupling constants in the range and manner defined 
above. 

The left panel of Fig. [2] shows the dependence of the 
EoS of hypernuclear matter on variations of the cou- 
pling constant in the range 0.26 < XY,Y.a < 0.66 at 
fixed XAh<y = 0.58 (a value that corresponds to the NSC 
model), where xhhu = guHa / gNNa is the ratio of the 
hyperon-cr and the nucleon-cr coupling constants. The 
resulting EoS cover the shaded area, which is bound by 
the two EoS corresponding to the limiting values of the 
parameter gsSff! in addition we show the special case 
where the ^seo- parameter is fitted to reproduce the em- 
pirical value of the S-hyperon potential in nuclear matter 
at saturation [/s = 30 MeV according to 



400 




gNNcr 



gwNL, 



[UJ 



(31) 



FIG. 1: (Color online) Equations of state of nuclear (dashed, 
blue online) and hypernuclear (solid, red online) matter. The 
coupling constants for nuclear matter correspond to density- 
dependent (DD) ME2 parametrization of Ref. [s^. The cou- 
pling constants of hyperons are related to the nucleonic ones 
according to the SU(6) symmetric quark model. The density 
is normalized to the nuclear saturation density po = 0.152 
fm^"' in the DD ME2 parametrization. 



where (...) refers to a mean- field value of the field. The 
right panel of Fig. [5] shows the dependence of the EoS of 
hypernuclear matter on variations of the gAAa coupling 
constant in the range 0.52 < XA\a 0.66 at fixed x^iscr = 

0. 448 (as implied by the NSC model). The set of resulting 
EoS covers the shaded area; we also show explicitly the 
limiting cases as well as the special case, where the (?aA(t 
coupling is fitted to reproduce the A-hyperon potential in 
nuclear matter at saturation Ua = —28 MeV according 
to an equation analogous to Eq. (pij) . As expected, the 
reduction of the hyperon couplings stiffens the EoS. The 
stiffest hypernuclear EoS in this set corresponds to the 
values of couplings xsso- — 0.448 and XAAa = 0.52. 

Fig. [3] shows the particle fractions of fermions, defined 
as the ratio of particle number densities to the baryon 
number density for the limiting cases of the SU(6) quark 
model couplings and the stiffest hypernuclear EoS. There 
are a number of features that are common to the both 
EoS. The appearance of S~ hyperons, which compensate 
the positive change of the nuclear matter and excess of 
negative isospin of the neutrons, is favored at low density. 
They contribute to the suppression and eventual extinc- 
tion of electron and muon populations at large densities, 

1. e., the matter becomes deleptonized. Once the increas- 
ing with density population of hyperons reaches the 
population of protons, the charge neutrality prevents fur- 
ther increase of S~ hyperons. Instead, an excess of neu- 
tral A hyperons builds up and dominates the hyperon 
population at asymptotically large densities. In the case 
of softer EoS (left panel of Fig. the dclcptonization 
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FIG. 2: (Color online) Equations of state of hypernuclear 
matter for a range of values of hyperon-cr meson couplings 
defined in terms of xhhct ~ gHH<y/gNNcr, H £ A,Y:. The nu- 
cleonic EoS (dotted line, magenta online) and hyperonic EoS 
with SU(6) quark model couplings (dot-dashed line, green 
online) are shown as a reference. The nucleonic coupling 
constants correspond to the DD-ME2 parametrization [s^ : 
the hyperon-vector meson couplings are fixed as explained 
in the text. Left panel: we assume xaact ~ 0.58, as in 
the NSC potential model, and a range 0.26 < xsect < 0.66 
which generates the shaded area. Right panel: we assume 
xsstT ~ 0.448, as in the NSC potential model, and a range 
0.26 < XAAa < 0.66. The cases xeect = 0.46 (left panel) 
and XAAtT ~ 0.59 (right panel), shown by dash-double-dotted 
(black online) lines, fit the depth of the potentials of the E~ 
and A hyperons in nuclear matter at saturation. 



and onset of hyperons like and cascades EP'~ oc- 

curs at densities lower than in the case of the stiff EoS 
(right panel of Fig. [31). This is the consequence of the fact 
that the large hyperon-meson couplings favor the forma- 
tion of hyperons, which in turn soften the EoS. Thus, 
the SU(6) quark model, with its large meson-hyperon 
couplings favors the onset of hyperons at low densities, 
which produces a soft EoS. Conversely, the stiffest EoS, 
which corresponds to small values of the hyperon-meson 
couplings, disfavors the onset of hyperons. 



B. Finite temperatures 

Compact stars become transparent to neutrinos after 
about a minute of formation in a supernova explosion. 
Prior to the neutrino-transparency era the star is hot, 
with temperatures of order of several tens of MeV and 
neutrinos are in equilibrium with matter. The neutrino 
thermal distribution is characterized by the neutrino 
chemical potential /i^ and neutrino fraction Fj, = u^/ub, 
i.e., the ratio of the neutrino number density to the 
baryon number density. The neutrino fraction is com- 
monly parameterized by specifying the Icpton fraction for 
each flavor = Yi + Yi^, where Y/ is the fraction of the 



FIG. 3: (Color online) Particle fractions in hypernuclear 
matter at T = 0. Left panel: hyperon-scalar meson cou- 
plings are fixed as in the SU(6) symmetric quark model. Right 
panel: a stiff hypernuclear EoS from our parameter space with 

T^^- = n ddS and » - = n .(^9 
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FIG. 4: (Color online) Equation of state of hypernuclear 
matter at T = 50 MeV with and without neutrinos. In the 
left panel the hyperon-scalar meson coupling constants are 
fixed by the SU(6) quark model, whereas in the right panel 
(/SECT = 0.448(;jvjvct and gAA<j ~ 0.52gNNrT, which corresponds 
to the stiffest cold hypernuclear EoS considered. In each panel 
the dash-dotted line (green online) corresponds to the case 
without neutrinos. The shaded region shows the range of 
variation of each EoS when the lepton fraction is changed 
between the limiting cases Yl —0.1 (dashed, red online) and 
Yl = 0.4 (solid, blue online). 



charged leptons of given flavor. The relevant to proto- 
neutron stars range of lepton fraction is 0.1 < < 0.4. 

Figure 3] shows the EoS of hot hypernuclear matter 
in the cases when there is a thermal population of neu- 
trinos with Q-l < Yl < 0.4 (shaded area) and in the 
absence of neutrinos, i.e., Yi^ ~ and Yi determined 
from /^-equilibrium with = 0, as in cold hypernuclear 
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FIG. 5: (Color online) Particle fractions in hypernuclear mat- 
ter at r = 50 MeV for a soft hypernuclear EoS with couplings 
fixed by SU(6) symmetric quark model; — in the left 
panel and Yl = 0.4 in the right panel. 



matter. The influence of the neutrino population was 
considered for two parametrizations: one corresponds to 
the parameters implied by the SU(6) quark model and a 
soft EoS and another corresponds to xsso- = 0.448 and 
XAAa = 0.52 and a stiff EoS (see Fig. [D for the T = 
counterparts). It can be seen that the EoS is stiff when 
neutrinos are present {Y^ ^ 0) and the larger is the frac- 
tion of neutrinos, i.e. Yl, the stiffer is the EoS. The 
stiffening of the EoS can be attributed to the fact that 
the thermal population of neutrinos adds its contribution 
to the pressure of matter. However, neutrinos stiffen the 
EoS also indirectly by changing the composition of mat- 
ter. As seen in Figs. [S] and [5] the presence of neutrinos 
has a strong effect on charged leptons: the deleptoniza- 
tion effect observed in neutrino-less matter is reversed 
when neutrinos are present. As a consequence the E~ 
hyperons are disfavored in neutrino-rich matter. In fact 
the sequence of appearance of charged S hyperons is re- 
versed (cf. the left and right panels of Fig. [S|) . Because 
neutrinos suppress the hyperon (predominantly A) frac- 
tions, the EoS is closer to the nucleonic one, which as we 
have seen, is always stiffer than the hypernuclear EoS. 
By comparing Figs. [S] and [5] one can access how various 
choices of the hyperon-scalar-meson couplings influence 
the composition of matter. The general trend, observed 
at T = (Fig. [3]), persists at non-zero temperatures, i.e., 
the thresholds for the onset of the various hyperons shift 
to higher densities for smaller values of the couplings of 
hyperons to the scalar mesons. The models with high 
values of these couplings are softer and more hyperon- 
rich as the ones with smaller values of these couplings. 



C. Polytropic fits to the EoS 

Numerical relativity and astrophysics problems fre- 
quently require a simple parametrization of the EoS of 
matter at zero temperature. The polytropic form of the 
EoS has been one of the common choices and, more re- 
cently, piecewise-polytropic EoS were used to construct 
accurate representations of microscopic EoS [6^. We 
shall use below the following representation 

4 

P = Y,K,p''^e{p~a,po)e{hpQ - p), (32) 

i=l 

where F^ is the polytropic index, Ki is a dimensionful 
constant, po is the saturation density. These parame- 
ters are used to fit the EoS within a density interval, 
whose lower and upper boundaries are specified by the 
parameters and bi. The six different EoS that were 
parametrized using Eq. (j32p include a nuclear EoS, a 
hypernuclear EoS based on the SU(6) symmetric quark 
model, and four additional EoS with x'sso- and x'aao- from 
the range discussed above. Table [Tl] lists the fitting pa- 
rameters. The polytropic representation above should 
be matched to an EoS of matter in the crust of a neu- 
tron star at the density Q.5po (assumed to be the largest 
density in the crust). To optimize the fit, the density seg- 
ments were chosen to optimally reflect the changes in the 
EoS; for example one knot, in the case of hypernuclear 
EoS, is chosen to be the density of the onset of hyper- 
ons. The change of the polytropic index with density is 
consistent with the expectation that as matter becomes 
more relativistic, the polytropic index decreases. 
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FIG. 6: (Color online) Same as in Fig. O but for a stiff 
hypernuclear EoS with couplings a-EEo- = 0.448 and a-AA^ = 
0.52. 



TABLE II: Piecewise-polytropic parametrizations of various 
EoS according to Eq. H32|) . Set A corresponds to a pure nu- 
cleonic EoS. The hyperonic EoS with different values xsi^a 
and XAAa are arranged as follows: set B: xaact ~ 0.52 and 
a;EECT = 0.448, set C: XAAa = 0.66 and xeeo- = 0.448, set 
D: XAAtT = 0.58 and x-sstr = 0.26, set E: XAAa = 0.58 and 
x-EStT = 0.66. Set F refers to the hyperonic EOS with xhhct 
fixed to the values of the SU(6) symmetric quark model. The 
dimensionful constant Ki is given in units of MeV ira^+^^- . 
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FIG. 7: (Color online) Dependence of the gravitational mass 
of compact hypernuclear stars on central density at zero tem- 
perature. The solid (blue online) and dashed (red online) 
lines show the limiting cases of parameter space as indicated 
in the panels. The dash-dotted (green online) line shows the 
observational lower limit on the maximum mass 1.97M0. 



D. Stellar structure 

The spherically symmetric solutions of Einstein's equa- 
tions for sclf-gravitating fluids are given by the well- 
known Tolman-Oppenheimer-Volkoff equations, which 
can be integrated for any given EoS. It is convenient 
to parametrize the equilibrium sequences of non-rotating 
configurations at zero temperature by their central den- 
sity. A universal feature of these solutions is the exis- 
tence of a maximum mass for any EoS, i.e., as the central 
density is increased a sequence reaches the configuration 
with the maximum mass and the stars with larger cen- 
tral densities are unstable towards gravitational collapse. 
A condition of stability for a sequence of configurations 
is dM/dpc > 0, i.e., the mass should be an increasing 
function of the central density. Alternatively, the sta- 
bility analysis of the lowest-order harmonics of pulsation 
modes (e.g. the fundamental radial pulsations) allows 
to access the stability of a configuration, for these are 
damped for stable stars and increase exponentially for 
unstable stars. 

The range of considered hypernuclear EoS trans- 
lates into the band of stable configurations shown in 
Fig. [T] The configuration containing hypernuclear mat- 
ter branch off from the purely nuclear ones as the cen- 
tral density is increased beyond the threshold density 
for appearance of hyperons in matter. The gravitational 
masses of hypernuclear stars increase with the density in- 
dicating a stable branch of these objects. The maximum 



masses are < 2.25Mq and are reached for central densi- 
ties of order 7po. Most of the sequences generated by the 
parameter space of the couplings considered are compat- 
ible with the observational bound M/Mq > 1.97. There 
is also room left for larger mass stars to allow for statisti- 
cal distribution of the masses of neutron stars beyond this 
limit. Note that IAMq canonical mass stars would be 
pure nucleonic for the stiffer subclass of EoS considered, 
whereas they would contain hypernuclear matter for the 
softer subclass; however for all the EoS considered, the 
stars with M > 1.5Mq contain hypernuclear matter. 

In Fig. [8] we show the mass-radius relationship for the 
hypernuclear sequences. The configurations with masses 
close to the maximum mass M ~ 2Mq have radii of 
i? ~ 12 km, whereas the canonical mass stars have radii 
of order 14 km. The figure also shows the bound, which 
predicts that PSR J0437-4715, being a Af = 1.76Mq neu- 
tron star, has a radius R > 12.5 km at 2(t level. Clearly, 
the hypernuclear stars are consistent with this observa- 
tion. Finally we note that it is not excluded that the 
low-mass neutron stars with M ~ 1.2A/0 may already 
contain hypernuclear matter. 



VI. CONCLUSIONS 

Despite of decades long theoretical research on the hy- 
pernuclear matter, the appearance of hyperons in com- 
pact stars remains an open issue. While the recent as- 
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FIG. 8: (Color online) The mass-radius relations for compact 
hypernuclear stars at zero temperature. The labeling and pa- 
rameter space is as in Fig. [T] The arrow shows the mass-radius 
constraint of Ref. [sj] at 2a level, which is M = 1.76Mq and 
R > 12.5 km. 

trophysical measurements exclude a significant fraction 
of soft EoS, the hyperonization of dense nuclear matter 
remains a serious possibility. Our present study confirms 
this within a relativistic density functional approach to 
nuclear matter, where we investigated the impact of vari- 
ation of the hyperon-scalar meson couplings on the EoS 
of hypernuclear matter. The range of found EoS is suf- 
ficiently stiff to produce heavy enough compact stars 
(M < 2.25M0). The radii of our sequences are located 
in the range of 12 < i? < 14 km. Piecewise-polytropic 
fits for six representative EoS are provided, which span 
the complete range of EoS from our parameter study. 

The parameter space of couplings of hyperons to scalar 
mesons was explored holding density dependent nucle- 
onic couplings fixed to their values suggested by the 
DD-ME2 parametrization of the nuclear density func- 



tional [56[. To allow for hyperonization in massive stars 
a requirement is to have small ratios of the hypernuclear 
to nuclear couplings, in particular, hyperons should be 
coupled to scalar mesons weaker than predicted by the 
SU(6) quark model. 

Our study of the EoS and composition of matter at 
non-zero temperature and including thermal cnsamble 
of neutrinos (present in a compact star during the first 
minute after birth) confirms that the neutrinos stiffen 
the high-density EoS and impact the charge lepton con- 
tent of hypernuclear matter. Instead of deleptonization 
with increasing density, seen in neutrino-less matter, the 
abundances of charged leptons remain constant, which 
has the consequence that the thresholds for appearance 
of charged E's is reversed. 

The hyperonic abundances found in our EoS are 
broadly consistent with the predictions of other (both 
early and recent) models. Indeed, these are determined 
by the hyperon mass spectrum and conditions of charge 
neutrality and /3-equilibrium. By varying the scalar me- 
son couplings, we showed that the stiffer the EoS the 
larger are the densities at which any type of hyperon 
appears. For example, substantial amounts of cascades 
appear already at densities of the order (3 — 4)po for a soft 
EoS and only trace amounts for a stiff EoS. Furthermore, 
at non-zero temperatures the thresholds for the onset of 
hyperons are located at densities lower than those at zero 
temperature. It is evident that during the early cooling 
stage of a neutron star (the time period when the tem- 
perature drops from tens of MeV to a MeV and matter 
becomes transparent to neutrinos) a rearrangement of 
particle content of matter must take place. 
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